The Three-Nucleon System at Next-To-Next-To-Leading Order 
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We calculate higher order corrections for the three-nucleon system up to next-to-next-to-leading 
within an effective field theory with contact interactions alone. We employ a subtraction formalism 
previously developed and for which it has been shown that no new three-body force counterterm 
is needed for complete renormalization up to this order. We give results for the neutron-deuteron 
phaseshifts and the triton binding energy. Our results are in very good agreement with experimental 
results and calculations using realistic nucleon-nucleon potentials. 
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I. INTRODUCTION 

The effective field theory (EFT) with contact interac- 
tions alone allows for a systematic calculation of low- 
energy two-body observables in terms of the effective 
range parameters 0] . With the correct power counting it 
can be applied to systems with large scattering length a 
and is an expansion in R/a, where R denotes the range 
of the underlying interaction. 

When applied to the three-body system with large scat- 
tering length, a three-body force [3| ~ or equivalently one 
subtraction - is needed to obtain cutoff independent re- 
sults at leading order (LO). Although this fact seems to 
limit the predictive power of the EFT, it has been very 
successful in describing various atomic and nuclear low- 
energy observables and has been used as a tool to un- 
derstand universal properties of few-body systems with 
large scattering length [1, 0, Q • 

To achieve high precision in three-body calculations 
higher order corrections have to be included. It is clear 
how to do this for the two-body system, however, in the 
three-body system it is not a priori obvious at which 
order an additional three-body counterterm has to be in- 
cluded. While Hammer and Mehen showed that no 
additional three-body datum is needed for renormaliza- 
tion at next-to-leading order (NLO) and later works have 
come to the same conclusion, different conclusion have 
been reached for different renormalization methods at 
next-to-next-to-leading order (NNLO). While Bedaque et 
al. and Barford and Birse found that an additional 
energy-dependent three-body counterterm is needed for 
renormalization if an explicit three-body force is used and 
the cutoff A is kept at ~ 1/R , we showed recently 
that within a subtraction formalism previously developed 
fldi [m , no additional three-body counterterm has to be 
introduced at NLO and NNLO for a consistent renormal- 
ization of observables if A ^ 1 /R. In this work, we com- 
puted also effective range corrections for the *He trimer 
system up to NNLO and achieved very good agreement 
with a previous calculation using a realistic atom-atom 



potential. We showed that observables are cutoff inde- 
pendent for sufficiently large cutoffs and that effective 
range corrections to observables scale as expected. In 
this paper we will extend this work to the three-nucleon 
system and discuss the corresponding results. As the 
power counting in the neutron-deuteron quartet channel 
is well understood jl^ , we will focus here on the neutron- 
deuteron doublet channel. 

In the following we will explain briefly the subtraction 
scheme which is used in this work to rcnormalizc three- 
body observables, give results for the triton binding en- 
ergy and scattering phaseshifts and discuss corrections up 
to NNLO to a universal correlation between the nucleon- 
deuteron scattering length and the triton binding energy, 
known as the Phillips line. We will end this paper with 
a short conclusion and an outlook. 



II. THEORY 

At sufficiently low energies few-body systems interact- 
ing through short-range interactions can be described 
with an EFT built up from contact interactions alone. 
Employing an auxiliary field formalism ^ the neutron- 
deuteron system in the singlet channel can be described 
by the Lagrangian |^ 
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where N represents the nucleon field and ti{sj) are the 
dinuclcon fields for the 5*1 (^6*0) channels with the cor- 
responding quantum numbers, respectively. The dots in- 
dicate more terms with more fields/derivatives. 
The full two-body propagator r is the result of dressing 
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the bare dimicleon propagator by nucleon loops to all 
orders 



(2) 



where the index a — s,t indicates either the singlet or 
tiplet two-nucleon scattering channel. In this form the 
the two-body propagator has poles at energies outside 
the validity region of the EFT. Therefore the propagator 
cannot be used within the three-body integral equation 
for cutoffs A > 1/r without employing additional tech- 
niques to subtract these unphysical poles. Instead of us- 
ing the propagator in the form above we will expand it 
up to a given order in R/a 



ri-\E)^'-^. 
^ ' E + Bd 
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with Bd — jf/m the deuteron binding energy and for 
n < 3 S'(") is defined as 

^^"'(^)-^^(D + (4) 

The set of integral equations for nucleon-deuteron scat- 
tering generated by this EFT (neglecting the three-body 
force for the beginning) is given by 0, [llj 

Ki:\q,q';E) = Zu{q,q';E) 
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where m denotes the nucleon, n the order of the cal- 
culation (for the following definitions it will always be 
assumed that n < 3) and Zafs the Born amplitude 
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with the iso-spin matrix A given by 
1/1 -3 
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The set of integral equations in Eq.Q is strongly cutoff 
dependent and a three-body force has to be introduced or 
equivalently a subtraction has to be performed to render 
observables cutoff independent. At threshold the integral 
equations are renormalized by noting that 



K[:\Q,Q--Bd) = 



and subtracting this known quantity from Eq.Q 
K['l\q,G--Bd) - Ki'^\0,0;~Bd) , 
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After rewriting the resulting set of integral equations 
the half-off-shell threshold amplitude takes the following 
form [IIIOI 

Ki:\q,0;-Bd) = K[-\Q,Q--Bd) 

+ [Zu{qA,~Bd)~ Zu{%0;~Bd)] 
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(10) 

The amplitudes Kj:f\q,0; —Bd) and isr^"''(g,0; —Bd) are 
fully renormalized after the subtraction is performed. An 
essential point in obtaining the amplitudes at any energy 
and momentum is demanding that 



Ki:^\q,0;-Bd) = K'^;"{0,q;-Bd) , 
Ki:\q,0;-Bd) = Ki'^\0,q;-Bd) . (11) 
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Using resolvent identities, subtracted integral equations 
for any energy can be derived. For further details, the 
reader is advised to turn to 0,^3- 

The two-body parameters used throughout the rest of 
this work are given by 



7r' = 4.317 fm , 
77^ = -25.04 fm 



n = 1.764 fm , 
r, = 2.73 fm . 



(12) 



We will use the result of a recent neutron-deuteron scat- 
tering length measurement with — 0.645 ±0.005 fm as 
our three-body input [l^. 





parameters 


as [fm] 


Bi [Mev] 


LO 




0.645 


8.08 


NLO 


a3,'yt,as,rt,rt 


0.645 


8.19 


NNLO 




0.645 


8.54 


EXP 




0.65 


8.48 



TABLE I: EFT predictions for the triton binding energy up to 
NNLO using the neutron-deuteron scattering length as three- 
body input. The second column indicates which parameters 
were used at what order of the calculation. Energies and 
lengths are given in MeV and fm, respectively. 



III. RESULTS 

The three-nucleon system has been previously consid- 
ered within the EFT with contact interactions alone. 
Leading order results for phaseshift and the correla- 
tion between triton binding energy and neutron-deuteron 
scattering length were obtained in Q . Bedaquc et al. 'T'I 
calculated higher order corrections up to NNLO intro- 
ducing an additional energy-dependent three-body force 
at this order. In this section we will demonstrate by 
explicit calculation that no additional three-body force 
is needed for renormalization and considerable improve- 
ment in results is achieved at NNLO. We have computed 
the triton binding energy at LO, NLO and NNLO us- 
ing the neutron-deuteron scattering length as the three- 
body input parameter Our results are displayed 
m table n When using the experimental value for the 
neutron-deuteron singlet scattering length 03 = 0.645 fm 
as three-body input we obtain a triton binding energy of 
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FIG. 1: (Color online) Phase shifts for neutron-deuteron scat- 
tering below the deuteron breakup at LO (dashed-dotted 
line), NLO (dashed hue) and NNLO (soUd line). The trian- 
gle is the result of the scattering length measurement of [l^ . 
The circles are the results of the van Oers-Seagrave phaseshift 
analysis p3 . and the squares denote a phaseshift calculation 
using a realistic nucleon-nucleon potential 



= 8.48 MeV at NNLO which has to be compared to 
the experimental value I^'^^ = 8.54 MeV. The results 
show a significant improvement with each order and are 
sufficiently close at NNLO to the experimental value to 
agree with a projected error at this order of approxi- 
mately (rt^t)^ ~ 3%. At these orders the error caused 
by the uncertainty in the effective range parameters is 
smaller than the error caused by N'^LO corrections. 
Using the set of integral equations in Eq. Q with the ex- 
panded two-body propagator in Eq.© leads at LO, NLO 
and NNLO to the phaseshifts shown in Fig^ For com- 
parison we show in the same figure the results of a forty 
year old phaseshift analysis Il4| and a theoretical calcu- 
lation using a realistic nucleon-nucleon potential |l5j . At 
higher order our results seem to describe the experimen- 
tal data better but considering the age of the analysis 
and the fact that no errors are given for these data, the 
theoretical calculation by Kievsky et al. should be con- 
sidered as the true benchmark test for our calculation. 
At NLO our results already lie significantly closer to this 
calculation and nearly perfect agreement is achieved at 
NNLO. It should be noted that our results at LO and 
NLO order agree with previous EFT calculations results 
given in FRlll| . We also achieve qualitative agreement 
at NNLO with Ref.0, however, without employing an 
additional three-body counterterm. 

A further way to illustrate the improvement in our re- 
sults is to consider the Phillips line. The Phillips line 
is a universal feature of three-body systems with a large 
two-body scattering length and arises as a nearly linear 
correlation between the 1-1-2 scattering length and the 
three-body binding energy. In Fig|21 we display our re- 
sults for LO, NLO and NNLO and display also the exper- 
imental value. It is interesting to note that in contradis- 
tinction to the ''He trimer system the Phillips line has 
not converged to a definite result yet 9j. The obvious 
reason for this is the rather large expansion coefficient 
which is roughly jti^t '^-^ 1/3 while jr ~ 0.1 in the ^He 
system. The power counting derived in jQi] is valid for 
cutoffs much larger than 1/r. Therefore, our results have 
been evaluated at corresponding cutoffs and in fact our 
results are fully converged with respect to A. The re- 
sults given in this section are also numerically converged 
up to the digits displayed. Further, we also analyzed 
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FIG. 2: (Color online) The Phillips hne for leading (dot- 
dashed line), next-to- leading (dashed line) and next-to-next- 
to-leading (solid line) order. The cross denotes the experi- 
mental value. 

the convergence of three-body observables when using a 
different parametrisation for hi gher order two-body cor- 
rections called Z-matching [T^ Il7j and have found the 
same results up to the level of accuracy expected in an 
NNLO calculation. 



results for the binding energies and scattering phaseshifts 
simultaneously. Our value for the triton binding energy 
agrees very well with the experimental value and the re- 
sults for the neutron-deuteron phaseshifts seem to be 
in nearly perfect agreement with a calculation using a 
realistic nucleon-nucleon potential. The results for the 
scattering phaseshifts a gree qualitatively with previous 
calculations at NLO [ITTH and NNLO 7}, ahhough a 
second three-body counterterm was included in 0. In 
particular, these results seem to agree very well with an 
expected error of (7r)3 ^ 3% for an NNLO calculation. 
We have therefore presented further numerical evidence 
which supports the claim in Ref. Q that calculations can 
be performed with exactly on three-body counterterm up 
to NNLO in the EFT with contact interactions alone. 
Our analysis presented in ^ indicates that at N^^LO an 
additional three-body input is needed for renormalization 
of observables. Therefore, NNLO can be also considered 
as the last order at which a prediction can be made for 
the Phillips line as it is a correlation gouverned by one 
three-body parameter. 

Further possible applications of the subtraction for- 
malism include the calculation of scattering observables 
above the breakup threshold and the coupling of exter- 
nal currents to the three-nucleon system, includiiig the 
electromagnetic form factor beyond leading order |l8| . 



IV. SUMMARY 

111 this paper we have computed three-nucleon observ- 
ables with the EFT with contact interactions alone up 
to NNLO using a subtraction scheme for renormaliza- 
tion. We have shown that the results for observables 
improve significantly at NNLO without performing an 
additional subtraction, i,e. without the introduction of a 
further energy dependent three-body counterterm. This 
improvement is in particular obvious if one considers the 
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